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ABSTRACT 

The numerical frontogenes is model of Williams [1974] is 
converted to Jin p vertical spacing and modified to include 
a higher upper boundary, tropopause and stratosphere, mois- 
ture and variable vertical turbulent diffusion coefficients 
of momentum and heat. The numerical solutions show more 
realistic quasi-steady fronts forming within 1-2 days. 

These solutions are examined and compared over a range of 
various boundary value coefficients and various ranges of 
temperature . 
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I. INTRODUCTION 



Hoskins and Bretherton [1972] have shown analytically, 
and Williams [1967, 1972] has shown numerically that dis- 
continuous fronts can form within a finite period of time 
if no turbulent diffusion is present. These studies sug- 
gest that a discontinuity will form within 24 to 36 hours 
when reasonable initial conditions are used. If turbulent 
diffusion is present, it can be expected that a balance will 
be achieved between the f rontogenetic advections and turbu- 
lent diffusions of momentum, heat, and moisture. The front 
should remain in this state of quasi -balance as long as the 
large scale deformation field causes frontogenetic advec- 
tions around the front. 

The numerical frontogenesis model of Williams [1974] 
shows the effects of horizontal and vertical turbulent dif- 
fusions of momentum and heat in the formation of these 
quasi-steady fronts. This model which includes an Ekman 
boundary layer will hereafter be referred to as W74 and the 
model as it existed in Williams [1972] will be referred to 
as W72. Cornelius [1974] added moisture with its subsequent 
condensation and release of latent heat to W74. This in- 
vestigation raises the upper boundary, includes a tropopause 
and stratosphere and appends variable vertical turbulent 
momentum and heat diffusion coefficients to the work of 
Cornelius [1974] . 
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The purpose of this study is to obtain and examine frontal 
solutions throughout a four-day period. The physical model 
is essentially the same as one of the models treated by 
Hoskins [1971] and Hoskins and Bretherton [1972] . The 
basic model, as it now exists, numerically describes frcnto- 
genesis which is forced by a nondivergent horizontal wind 
field which contains stretching deformation. This deforma- 
tion wind field is constant in time and independent of height, 
except in the Ekman boundary layer. The moist hydrostatic 
primitive equations with diffusion are used in which the 
time dependent quantities are functions of y and p only. 
Therefore, output from the model is in the form of vertical 
cross-sections in the Y-P plane. 

The model of W74 has been transformed to pressure coordi- 
nates to obtain more accurate values of pressure for the 
moisture equations of Cornelius [1974] . 

In Section II, the basic forecast equations are developed 
and a simplification which keeps the problem two-dimensional 
is discussed. The initial conditions are presented in Sec- 
tion III. Convective adjustment processes and latent heat 
release in the model are included in Section IV. The experi- 
ments undertaken are examined in Section V and conclusions 
are given in Section VI. 
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II . BASIC EQUATIONS AND BOUNDARY CONDITIONS 



The hydrostatic primitive equations with diffusion and 
moisture may be written as follows: 



3V 3 (wV) 

_ + V* (VV) + -Jp- + V<f> + f(kxv) = A m V 2 V 

3 

+ 3p (K M Fp ^ ' 

|| + v- (ev) ♦ j| ( M 0) - AgV 2 0 * j| (X 0 || ) 

+ Q A + Q M + H c > 

If + V -W - fj- ( M q) - A q V 2 q ♦ |p Of, If ) 
+ m a - M c • 

V-V + = 0, 

3p 

If = " L( p } > 



where 

L(p) 



= R r 2_ 



= c 



p 3i 






( 2 . 1 ) 



( 2 . 2 ) 



(2.3) 

(2.4) 

(2.5) 



The equations are in pressure coordinates with slight modi- 
fication due to the addition of the equation of conservation 
of water vapor (2.3) which is analogous to the thermodynamic 
equation (2.2). Sinks and sources of heat and moisture have 
been added to the thermodynamic equation (2.2) and the con- 
servation of water vapor equation (2.3). 
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The quantities, A^, Aq and A^ are the horizontal turbu- 
lent diffusion coefficients of momentum, heat and moisture, 
respectively, while K^, and are the corresponding ver- 
tical coefficients. These coefficients were assigned various 
reasonable values in the following experiments. Dry and 
moist convective adjustment processes are represented by 
the functions and while latent heating due to large 
scale condensation is represented by H c> The quantity 
is a sink or source of moisture due to moist convective 
adjustment in the moisture equation, while is a sink for 
moisture due to large scale condensation. 

Boundary conditions for the above equations are such 
that the domain has a flat lower boundary and a constant 
pressure surface where w = 0 as the upper boundary. The no 
slip condition is used at the lower boundary and the no 
stress condition is used at the upper boundary. The condi- 
tions are as follows: 

w(x,y,0,t) = w(x,y,p t ,t) = 0 

V(x,y,0,t) = 0 

sv 

K\| "gT)" ( x >y>p+->t) - o 

M dp Z (2.6) 

K e If O-y" 0 - 1 ) - K e If (x,y,P t ,t) - o 

K q Tp Cx.y .0 ,t) - K q |3. (x,v,p t ,t) = 0 

where p represents the pressure level at the top of the 
domain, and 0 corresponds to p^ at the bottom. 
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The approximate steady-state solutions to the hydrostatic 



primitive equations [(2.1) - (2.5)] are: 



V = U = D [x(l-e ^cos £) + ye ^ sin £]i_ 

- r - r 

y(l-e ^cos c) - xe sin C)j_, 



(2.7) 



to = 0 , 



<{> = $ = -D 2 (x 2 +y 2 )/2 - fDxy + F(p), 
0 = G(p), 

q = H(p), 



where 




'M 



K,, is evaluated at p , D and £ are 
M *o ’ 



constants and where F(p) is such that the hydrostatic equa- 
tion (2.5) is satisfied. The functions G(p) and H(p) are 
arbitrary functions of p since to = 0 . These solutions sat- 
isfy the boundary conditions (2.6) and allow for inclusion 
of the Ekman boundary layer. The relations (2,7) satisfy 
the steady-state equations if some small advection terms in 
the boundary layer are neglected as they usually are in 
Ekman theory. 

The equations and the model become two-dimensional as in 
W74 and Cornelius [1974] after a subdivision of dependent 
variables as follows: 



V = U(x,y,p) + u(y,p,t)i + v(y,p,t)j_ 
“ = w(y,p,t) 

6 = e(y,p,t) 

4> = <J>(x,y) + tt (y,p,t) 
q = q(y»p»t) 
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It is noted that all departures from the approximate 
steady deformation solutions (2.7) are assumed to be inde- 
pendent of x. If we substitute the relations (2.8) in the 
component of (2.1) we obtain: 

+ fy (uv) + |p (wu) + uD (1-e^cos c) - D[y (l-e'^cos c) 

- xe ^ sin c] + vDe ** sin r, + wD |^- [xe~^(l-cos 5 ) 

+ sin 5 ]- fv + A M 1^ + (K m ). (2 . 9) 



In this study we will neglect the terms in (2.9) where 
x appears and apply this and the other equations at x = 0. 

If the initial u, v, w, 0, tt and q fields are independent of 
x, the error from this approximation will grow slowly and 
will be confined to the boundary layer since the x dependent 
terms of (2.9) are zero outside the boundary layer. Since 
frontogenesis occurs very rapidly, it is expected that the 
development of x-variations in the dependent variables would 
have only a small effect on the resulting quasi-steady front. 

In any case these effects could not be observed in the 
atmosphere, because variations in the basic fields along 
atmospheric fronts would be much more important. 

When equation (2.9) is evaluated at x = 0 it becomes 



3u 3/" -v . 9 ,• -v r 3V . 3u 

9"t Fy ^ UV -* Tp ru 3y rV ITy 



v fy sin ^ 



V ITp sin c) = fv + A m gyy + jp (K m F p ) , (2 . 10) 

where V = -Dy and T = 1-e^ cos £. The j_ component of (2.1) 
applied at x = 0 takes the form 
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5T + ly (vv) + fp (MV) * r 57 (vV ^ ‘ u 57 e "' sin 5 



- t«V) |p (e' ? cos O - - |y - fu ♦ A m 8y 

tK ™ 1 
^ K M 3p ' * 



3 2 u 3 

“ r 3p 



( 2 . 11 ) 



Similarly the following equation for the departure of poten- 
tial temperature can be derived from (2.2): 



li + 1 

3 1 



37 Cve) ♦ |p M) * rv §i 

+ + Q M + H c 



a 3 6 _ + 3 si, 3 6 

A e 3y2 3p ^ K 0 3p -* 



( 2 . 12 ) 



and the following equation for q from (2.3) 



ft + §7 * 3p 



(<oq) + TV I 5 - = A (K ^3. ) 

^ 3y q dy 2 3p v M 3p J 



+ m a - M c* 



(2.13) 



Equations (2.4) and (2.5) become 



and 




3 

3P 



= 0 , 



3 TT 

3p 



-L(p)6. 



(2.14) 



(2.15) 



The boundary conditions at the bottom and top of the domain, 
respectively, are 



U = V = 


= m = K q 


30 

3P 


K 

q 


= 

3P 


0 


II 

O 


v 9u . 
M 3P 


If 9v 
M 3p 


= 0) = 


K o 


3_e 

3p 


K 

q 


If - 0 p 



If we define the vertical average of a quantity 




(2.16) 
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1 



(2.17) 



< ( ) > = 



top 



P 

) v/ 

top 



)(dp) , 



integrate the hydrostatic equation (2.15) with respect to p 

and remove the vertical mean, it follows that 

p p 0 (2.18) 

7 T - < IT > = [ f L(p)0dp - < f L(p)0dp >] . 

Ptop Ptop 

Taking the vertical average of (2.14) and using the boundary 
(2.16) gives 



|y < v > = 0. (2.19) 

This equation states that the total mass flux in the y-direc- 
tion is independent of y. In W72 a symmetry argument was 
used to show that this flux must be zero. This argument 
does not hold strictly in this experiment, but it can be 
expected that at large distances from the frontal zone, the 
disturbance mass flux will vanish. Thus it is assumed that 



< v > = 0 . 



( 2 . 20 ) 



Equation (2.11) can be written 



& ♦ [w - <w>] ♦ |^( M v) * (rvv-<rv>v) 



[ue ^ sin z - <ue ^ sin z,> - V[w - 5 — (e cos ?) 

dp 



3y 



<w Tp cos 



fy (7T-<TT>) - f (U-<U>) + A M fy 



3 3V . L. , 3V , 

3p M 3p — ^ 3 p ^P=P 0 



( 2 . 21 ) 



The equations (2.10), (2.12), (2.13), (2.14), (2.18) and 
( 2 . 20 ) form a complete set which can be solved by a pure 
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N>|< 



inarching process. The finite difference equations conserve 
mean squares in the advection terms and are described by 
Williams [1967] . The Matsuno scheme is used everywhere 
except in the arrangement of variables and the finite dif- 
ference approximations are the same as those used by Williams 
[1967] . In order to close the problem, computational boun- 
daries must be introduced in y. Since the disturbance 
velocities should die out at a sufficient distance from the 
axis of dilation, then 

v(±y,p,t) = 0. (2.22) 

However, there is appreciable inflow across these computa- 
tional boundaries since V(±Y) = ± DY. The quantities u, 6 
and q which are advected across the boundaries must be spe- 
cified independent of the interior values if computational 
stability is to be maintained (Platzman, 1954). Thus, the 
following boundary conditions are used: 

u[±(Y + Ay/2), p,t] = u[±(Y + Ay/2),p,0] 

0[±(Y + Ay/2), p,t] = 6 [ ± (Y + Ay/2),p,0] 

(2 23) 

0[±(Y - Ay/2), p,t] = 0 [ ± (Y - Ay/2) ,p , 0 ] 

q[±(Y + Ay/2), p,t] = q [ ± (Y + Ay/2),p,0]. 

The computational boundaries y = ±Y are placed between grid- 
points so that the above conditions are actually applied at 
y = ± (Y + Ay/2). 

The boundary conditions were found to be satisfactory 
except that a solution separation developed near the boun- 
daries. This was controlled by using a Matsuno time dif- 
ferencing scheme for stability and because of its damping 
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characteristics. The vertical gridpoints were spaced pro- 

— s V 

portional to £n p using the relation P = P Q - p Q e , 
where s is Ap/H. Refer to Table II for forecast levels. 



21 



III. INITIAL CONDITIONS 



The initial conditions of Cornelius [1974] were modified 

by raising the upper boundary to approximately the 100 -mb 

level, and by simulating a "tropopause" centered at about 

the 250-mb level. The Cornelius [1974] initial conditions 

are developed as follows. The initial potential temperature 

field in pressure coordinates is given by 

c 9 0 " 

e(y>p>o) = 0 iCp o ) ex P [ -g — 9 y (i ■ ( ^— ) K ] 

- a( 2 /ir)arctan(sinh ay) (3.1) 

where 



a = 



f TT 



fn 



n/LCPn) (-3 0j/3p) Il/(R/p M ) (p/p M ) K (' 36 T / 9 P) 



(3.2) 



M J 



M' 



M 



In the model, a is given a constant median value by solving 
(3.2) at p = pj^ and using the static stability 36j/3p at 
this level. 

The initial u component of the velocity is found using 
a form of the hydrostatic equation (2.26) and the definition 
of geostrophic wind 



u 



1 3jf>_ 
f 3y ‘ 



A form of the thermal wind equation is found to be 



(3.3) 



3u _ 13 3<j> _ 1 R ( p 36 { , . n 

yp“"f373p-fp(Jr) 37- t3.4) 

We again assume a median pressure level where u(p^) = 0, 
and integrate 
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(3.5) 



r P 9 U a 
P J ^ dp 



R 39 p 1 ( P ^ k i 



to get 

u(y,p,0) 

where 



J2. ££ rr R_ ) K 

£ 3y 1 ‘- p Q } 



CfT )"] 

F o 



(3.6) 



K = R/c . 
p 

From equation (4.3), in Williams [1972], it is seen that 
3 0 2aa 



3y u sech (“y)- 

Therefore, the initial u field is given by 



(3.7) 



u(y,p,0) = - 2|« [( E_ . ( Em ) k ] sech (ay). (3.8) 

"o ^o 

The initial divergent wind is determined as in Williams 
[1972] by using the quas i -geostrophic equation in pressure 
coordinates. The stream function is defined as 



and 



v = dip/dp 



to = - 3i]j/3y 



(3.9) 



R 



3p 



M S K 



Nr 3y 



and equation (4.10) in Williams [1972] takes the form 

d z ip . T w 30 I , d 2 il> _ 2D T , 30 

3p2 L(p^)(- ) n 3v f2 ^(Pm) 3v » (3.10) 

recalling that 

l (Pm) = 



P M ^ 



M ' Pm ' p 



( ^ ) 
o 



The equation (3.10) is solved for ip from the initial temper- 
ature field (3.1) with the technique of Ogura and Charney 
[1962] and ip is then determined from (3.9). 
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These initial conditions do not satisfy the boundary 
conditions on u, v and 0 at the upper and lower boundaries. 
Hence a period of adjustment will be required to form the 
surface friction layer and some oscillations may be observed 
later. The fields of w and tt are obtained through the con- 
tinuity and hydrostatic equations (2.14) and (2.15). 

The initial field of specific humidity is calculated 
by first determining saturation vapor pressure e at each 
gridpoint through use of an integrated form of the Clausius- 
Clapeyron equation. According to Hess [1959] 



e = .611 
s 



exp ( 



R 



r L. _ 

1 T T 



k )) 



(3.11) 



V o 

where L is the latent heat of condensation, R is the gas 
constant for water vapor, and T q is a constant of integra- 
tion. With knowledge of the saturation vapor pressure, 
the saturation specific humidity is obtained from the ap- 
proximate relation of Haltiner and Martin [1957] 

.622 e 

q s = . (3.12) 



Knowledge of potential temperature at t = 0 from equation 
(3.1), makes it possible to calculate q g at each point using 
equations (3.11) and (3.12). Defining q to be a given 
fraction of q g allows the testing of atmospheres of different 
constant initial relative humidity. The initial conditions 
of Cornelius [1974] were modified as follows. A tropopause 
was simulated centered at about the 250 -mb level. This was 
accomplished using a reference temperature at the midpoint 
of the horizontal scale and calculating temperatures such 



24 



that a tropopause with its expected variation in height was 
determined and a zero vertical temperature gradient was 
achieved in the lower stratosphere. The tropopause height 
varied from 325 to 110 mb over the horizontal domain for 
0j (p ) = 290 K. The tropopause was simulated as follows. 

A reference temperature was computed using equation (3.1) 
as 

T re f = 215 + a ( | ) arctan(sinh ay) ( - ° -- Pr . ef ) K (3.13) 

^o 

When the temperature in each column reached this reference 
temperature, it was held constant at that value, as expected 
in the lower stratosphere. Temperature was then converted 
back to potential temperature. 
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IV. CONVECTIVE ADJUSTMENT, LATENT HEAT AND PRECIPITATION 



Two additional heating terms (Q m and H ) were included 
in the moist thermodynamic equation (2.2) necessitating 
further amendment to W74. Atmospheric heating by the dry 
convective adjustment, Q , was included in W74 to eliminate 

3 . 

any intense grid scale convection which developed due to the 
formation of any dry unstable lapse rates. This redistribu- 
tion of heating was accomplished using averaging techniques 
to ensure that the dry adiabatic lapse rate was not exceeded 
(3 0/ 3p < 0) . 

Lapse rates were tested in each vertical column of the 
model as in Cornelius [1974] . Upon encountering an unstable 
layer (30/3p > 0) , a cumulative average was taken and applied 
to the gridpoints within this layer until neutral stability 
(30/3p = 0) was attained. The model was checked simultaneously 
for instabilities developing below the convectively adjusted 
layer and potential temperatures at these levels were also 
cumulatively averaged and reassigned the new values. The 
final result of the technique is that the heating is redis- 
tributed to simulate physical processes of the atmosphere 
and mathematically that 30/3p< 0 throughout each column of 
the model. A completely analogous scheme was used for satu- 
rated parcels in the moist experiment, only using equivalent 
potential temperatures, 0g^, in place of potential temperature 
0 . 
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The addition of moisture required that large scale con- 
densation and the subsequent release of latent heat be given 
important consideration. Following Haltiner [1971], changes 
in moisture and temperature due to large scale condensation 
were determined as follows. 

Where relative humidity exceeds 1001, the adjustments in 
temperature 5T and specific humidity 6q are determined by 
the following equations: 



q + 6q = q s (T + 6T , p) (4.1) 

and 

c p 6T = - L c 6q. (4.2) 

These equations require that the final state of the air be 
exactly saturated and that excess moisture be condensed iso- 
barically, releasing the latent heat to the air. 

Equations (4.1) and (4.2) can be expressed in iterative 
form as 

9q 

Aq = q s (l - r) + ( -gf ) p AT, q = rq g (4.3) 

and 



c p AT = - L Aq 



(4.4) 



Eliminating AT from (4.3) and (4.4) gives 

q s (l ' r) 

Ac l = -T5q 



1 * ^ C It ) 



(4.5) 



At constant pressure the Claus ius -Clapeyron equation can be 
expressed as 



( 



3q s 

9T J p 



L q 
c n s 

PT ‘ 

V 



(4.6) 
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Substituting (4.6) into (4.5) gives 



q s Cl - r) 



(4.7) 




Aq may be now calculated from (4.7) as can the corresponding 
value of AT from (4.4). The first estimates of the adjusted 
temperature and specific humidity are, 



These values in turn are used to obtain improved approxima- 
tions. Iterations are continued until the values of T' and 
q* at each level match with sufficient accuracy 1001 rela- 
tive humidity. The changes are then applied to the potential 
temperature and specific humidity fields of the model ac- 
counting for the term of equation (2.2) and M c of equation 



The model is then investigated for unstable moist convec- 
tion and its associated heat of condensation. A convective 
adjustment scheme similar to the dry one is used on the satu- 
rated parcels so that moist static stability is always main- 
tained. A field of equivalent potential temperature is 
formulated according to Holton [1972] as 



and the moist convective adjustment is designed to ensure 
that the moist adiabatic lapse rate is not exceeded 



T' = T + AT and q' = q + Aq. 



(4.8) 



(2.3). 



0 



SE 



L q 

= 9 exp ( ) 

P 



(4.9) 



(80 SE /3p < 0). 
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The heat release, Q , as a result of moist convective 

^m 

adjustment is included through knowledge of the change of 
equivalent potential temperature during the adjustment. 

Taking the natural logarithm of both sides of equation 
(4.9) gives 



£n 0 SE 



Zn 0 



L q 
c n s 

c T 
P 



(4.10) 



A constant pressure surface may be assumed since the model 
uses pressure coordinates. Differentiating both sides of 
(4.10) and recalling the definition of potential temperature 



gives 



d0 

"0 



SE _ dT 
T 

SE 1 



dq L L q dT 

1 s c n 1 c 



c L s 



c p' r 



(4.11) 



Solving equation (4.6) for dq g , substituting into equation 

(4.11) and solving for dT leads to 

d0. 



dT = 



'SE 



©SE t^ 1 + R c^T 2 
v p 



^s L c 2 ^ 



(4.12) 



C P T 



where 

0 ). 



d0g E is the change in 0g E needed to adjust (90g E /3p = 



Since the atmosphere remains saturated throughout the 
process of moist convective adjustment, the Clausius -Clapeyron 
equation is again used to obtain the corresponding dq or 
actually dq g . Changes are again applied to the potential 
temperature and specific humidity fields of the model account- 
ing for the terms of equation (2.2) and M c of equation 
(2.3) . 
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A cumulative amount of precipitation is then found using 



the relation 



P 

in agreement 
condensation 
is expressed 




( 4 . 13 ) 



o 

with the changes in specific humidity due to 
and moist convective adjustment. Precipitation 
in inches of rainfall in the model. 
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V. NUMERICAL SOLUTIONS 



All numerical results to be shown use the following 
values for the constants: 



At = 540 sec 
Ap/H = .01978 

g = 9.81 m sec 



0 = 300 K 

o 



R = 461 joules 
Kg ' 1 K ' 1 



a = 12.56 K 



-2 



p = 100 cb 
*o 

Y = 1800 KM 
f = 10 4 sec 1 



p.. = 55 cb 
M 

Ay = 60 KM 
D = 10' 1 sec" 1 



0 T (p ) = 280,290,300 K c = 1003 joules 

P Kg ' 1 Kg 1 



R = 287 joules Kg 1 K 1 



k = R/c 



RH = 7 0 ! 



L c = 2.5xl0 6 joules 

v " 1 

K g 



ae 



3z 



1 = 5 K Km' 1 



ae T 

4P P M 



.53448 K cb' 



A m = A = 3xl0 4 m 2 sec 
m o 



0^ = 0, cb 2 sec 



- i 



A = 100000 m 2 sec" 1 

q 



C m B = .001409, .0007 cb 2 sec 



- 1 



(which corresponds to 0, 10, 5 m 2 sec 1 in height coordinates) 



. i 



'0B 



= 0, .0007, .0001761 cb 2 sec (which corresponds to 



- i 



0, 5, 1.25 m 2 sec * in height coordinates) 
where more than one value means that the constant was varied 
in different experiments. Whenever possible, the values are 
the same as those used in W74. 

The values of vertical diffusion of moisture C and heat 

m 

C 0 were varied in the experiments to compensate for density 
variation in the atmosphere. This is equivalent to a constant 
coefficient in the z coordinate system and is accomplished 
using the following formulae: 
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(5.1) 



C = C D (1 - p/p ) 2 
m mB ^ tf o J 

and 

c e “ C 6B (1 • P/Po) 2 (S- 2 ) 

in this: section numerical frontogenesis solutions with mois- 
ture, friction and a stratosphere are examined at t = 4 days 
and compared to the corresponding dry atmosphere. 

As can be seen from Table I, 11 experiments were performed. 
Experiments 1 and 2 contain only horizontal turbulent diffu- 
sions of momentum, heat and moisture while Experiments 3 
through 11 contain both horizontal and vertical turbulent 
diffusions of these quantities. 

In Table I, ®j(P 0 ) is the reference potential temperature 
(K) . The condensation interval is 5 time steps between suc- 
cessive applications of condensation and moist convective 
adjustment to the model in the moist experiments. Horizon- 
tal turbulent diffusion coefficients and initial boundary 
values of vertical turbulent diffusion coefficients (C m g, 

Cqb> Cq) are relative humidity (percent) and initial static 
stability (K/cb) . 

Experiments were performed with a reference potential 
temperature (0^^ Cp q ) ) > of 280K, 290K and 300K to simulate 
physically realistic situations. Changes to the dry atmo- 
sphere with different reference temperatures were almost 
negligible, but due to the exponential function in the 
Clausius.-Clapeyron equation important differences were ob- 
served in the moist atmosphere. 
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Experiments were also conducted using different boundary 
values of vertical turbulent diffusion coefficients of momentum 
and heat. 

The final experiment included a turbulent mixing coeffi- 
cient which reduced the values of and C Q g linearly to .1 

of their original values between approximately the 93- and 
84-cb levels in five equal steps. 

The adjusted values of vertical turbulent diffusion co- 
efficients of momentum C and heat C 0 were then varied as 

m 8 

described by equations (5.1) and (5.2). 

Primary emphasis on analyses of experimental results was 
on contour plots from the IBM 360/67 computer. These results 
obtained for all 11 experiments are analyses at t = 4 days 
of the fields of u, q, 0 and w. The entire vertical plane 
is shown from 100 to 10 cb in these plots. However, the 
outer portion of the domain between |y| = 1200 km and |y| = 

1800 km is not shown. 

Figures 1-3 show the initial fields of u, 8 and w at 
reference potential temperature of 290 K used in the follow- 
ing experiments. Contour plots for experiments 1 and 2 are 
shown as figures 4-10. Figures 4 and 5 show distributions 
of the rotational wind, u for corresponding dry and moist 
experiments, respectively. Features of u at all levels dif- 
fer quite significantly from dry to moist experiments. There 
is a marked increase in horizontal and vertical shear in the 
easterly flow in the moist experiment in the vicinity of the 
front near the surface. Horizontal and vertical wind shear 
is also more intense at low middle to middle levels in the 
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moist experiment. It is noteworthy that this property is 
dominant throughout the moist solutions. In the moist ex- 
periment the jet core increases in magnitude while shear 
decreases to zero at a lower level in the northern region 
near the tropopause. Figures 6 and 7 contrast values of 
6 throughout the domain. Noted is the warmer air in middle 
levels in the frontal zone. A slightly more intense gra- 
dient is seen near the surface with little change in tilt 
of the front. All this is attributed to the latent heat 
release in the moist experiment 2. 

Vertical motions for the experiments are investigated in 
Figures 8 and 9. Upward vertical motions occupy a much 
larger zone, extend much higher and occur closer to the 
frontal zone in the warm air mass of the moist experiment. 
Significant downward motion occurs between |Y| = 1200 and 
! Y | = 1800 km (not shown in the figures). 

Figure 10 shows the field of specific humidity for ex- 
periment number 2. Maximum precipitation is .37 inches. 
Figures 11 through 17 compare the dry and moist fields of 
u, 0 and w and show the q fields of experiments 3 and 4. 
Figures 11 and 12 show a small increase in the zone of 
easterly component in the lower levels of the moist experi- 
ment. Vertical and horizontal shear are also seen to in- 
crease as in experiments 1 and 2. Also noted is a much 
larger jet core in the moist experiment, with a slight 
northward shift. Noteworthy is a significant decrease in 
magnitude of the easterly component of u near the surface 
from that of experiments 1 and 2. Figures 13 and 14 exhibit 
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the same contrast in potential temperature as that of ex- 
periments 1 and 2. It is noted that the air at low levels 
and near the tropopause is colder in the moist experiment 
with vertical turbulent diffusion added (exp. 4). This is 
also true in the dry experiments. This is attributed to 
the increased amount of heat loss due to the vertical tur- 
bulent diffusion of momentum and heat offsetting the release 
of latent heat in the moist experiments. The air at mid 
levels in the frontal zone is warmer in the moist experi- 
ment as seen before. Figures 15 and 16 show a marked change 
in vertical motion between dry and moist experiments. A 
large zone of upward motion is seen near the frontal zone 

- 3 

with maximum values exceeding 3x10 mb/sec near the 85-cb 
level. Figure 17 shows the specific humidity of experiment 
4. Noteworthy is the dramatic shift to the south in specific 
humidity following the frontal position. 

Figures 19 through 38 compare the moist and dry fields 
of u, 0 and w and show the specific humidity of the moist 
experiments using boundary values of vertical turbulent dif- 
fusion of heat and momentum of .0007 cb 2 sec 1 (5 m 2 sec 1 

in height coordinates) at reference potential temperatures 
of 290, 280 and 300 K. 

Figures 18 and 19 contrast the u field of experiments 5 
and 6 at 0"j(P o ) = 290 K. Noteworthy is the increase in zone 
of maximum easterly component of the wind at low levels and 
the more intense horizontal and vertical shear in low levels 
over experiments 3 and 4. A slight northward shift and in- 
crease in size of the jet core is noted in the moist experiment. 
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Figures 20 and 21 show little difference in potential tem- 
perature fields, although the air in the frontal zone in 
middle levels is warmer in the moist experiment due latent 
heat release. There is a marked increase in tilt of the 
front over previous experiments, due to more normal n.orth- 
ward position of the front at the surface. Figures 22 and 
23 contrast the vertical motion of experiments 5 and 6 show- 
ing a maximum in upward motion near the 850-mb level con- 
sistent with previous experiments. Figure 24 shows the 
specific humidity field which differs slightly more in the 
same manner as experiment 4 differed from 2. Maximum pre- 
cipitation is .4 in. 

Figures 25 and 26 show more intense horizontal and ver- 
tical shear apparent in the easterly component of the moist 
experiment at 0 ^ Cp 0 ) = 280 K. The jet core is reduced in 
size and magnitude from experiments 5 and 6 and the hori- 
zontal shear decreases significantly near the tropopause. 

It should be noted here that the tropopause in experiments 
7 and 8 is lower than the tropopause in previous and subse- 
quent experiments due to the lower reference potential tem- 
perature used. Figures 27 and 28 show little difference 
between moist and dry fields of potential temperature except 
for the warmer air in mid levels in the frontal zone of the 
moist experiment due latent heat release. Noteworthy in 
these experiments 7 and 8 is the zone of low stability and 
apparent break in the tropopause in the northern part of the 
region. Figures 29 and 30 show a similar contrast between 
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moist and dry fields of vertical motion as, seen in previous 
experiments except that the contrast is not as dramatic as 
seen before. Figure 31 shows the q field is reduced con- 
siderably over previous experiments due to the lower refer- 
ence potential temperature used (280 K) . Maximum precipita- 
tion reduced to .15 inches. 

Figures 32 through 40 show a much higher tropopause 
due to the larger reference potential temperature (300 K) . 

A dramatic change in the u field is seen in Figures 32 and 
33 between dry and moist fields. Especially noteworthy 
is the very strong jet core and much more intense horizon- 
tal and vertical wind shear. Also note a large shift in 
the jet to the north which was not quite so obvious in pre- 
vious experiments. Figures 34 and 35 also show a marked 
contrast between dry and moist potential temperature fields. 
The air is much warmer in the upper levels near the tropo- 
pause due to the latent heat release. Thermal gradients 
along the frontal zone are also stronger in the mid and up- 
per levels. Strong heating in upper levels is appcirent. 

- 3 

Vertical motions (Figures 36 and 37) exceed 3x10 mb/sec 
in the warm moist air. The contrast is similar to that seen 
in previous experiments. The moisture field (Figure 38) is 
greatly intensified over previous experiments and specific 
humidities are in excess of 19 g/kg at low levels to the 
south. This is no doubt due to the larger quantities of 
latent heat release observed in this experiment, and the 
larger initial reference potential temperature. Maximum 
precipitation increased to .7] inches. 
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Figures 39 and 42 show experiment 11 which included the 
turbulent mixing coefficient into experiment 4. The only 
differences noted were a small pocket of lower specific 
humidity (Figure 42) at the top of the mixing level and a 
decreased zone of maximum upward motion in the warm moist 
air of the mixed experiment (Figure 41) . From the above 
we conclude that vertical turbulent mixing has no effect 
on mid and upper level deformation fields. 

Comparison of the results of the previous experiments 
with similar experiments of Cornelius [1974] revealed the 
following : 

1. The westerly wind field extended much lower while 
the easterly wind field was compressed below the 85 cb 
level. The zero wind level was also observed to be lower 
in the expanded atmosphere with tropopause included. This 
is true in both moist and dry experiments. The jet was also 
lower and more apparent. 

2. The pockets of warm air were apparent in the frontal 
zone as in Cornelius [1974] , but frontal tilt was somewhat 
less after 4 days in the experiments conducted in this study. 

3. Significant downward vertical motion occurred outside 
the zone observed in both studies although larger values are 
evident in the studies of Cornelius [1974] . 

4. Specific humidity fields were similar in both studies. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The numerical model of Williams [1974] is modified in 
this study to include moisture, a tropopause and lower stra- 
tosphere and variable vertical turbulent diffusion of heat 
and momentum. This leads to a more realistic situation and 
consequently more realistic solutions. Intensification of 
the baroclinic zone is seen to occur in upper levels, al- 
though no effect is apparent on the surface due to the in- 
clusion of moisture. The effects of vertical turbulent 
diffusion is apparent throughout the frontal zone through 
increased tilt, reduced heating and reduced values of wind- 
speed. The expanded boundary and inclusion of the tropo- 
pause along with a conversion to Jin p vertical spacing 
tended to shift all features to a lower level in the domain. 
Jet cores were generally larger and vertical motions in- 
creased in the warm moist air near the frontal zone due to 
the heating. The moist models illustrate an exponential 
thermal dependence in the various experiments at the reference 
potential temperatures of 280, 290 and 300 K. 

In future studies using this model, attention should be 
directed toward the use of still more realistic boundaries, 
such as a higher upper boundary. A translational motion 
should be added to the system so that warm and cold frontal 
effects could be investigated. Finally, an improved convec- 
tive scheme such as that of Arakawa and Schubert [1973] 
could be applied to the model. 
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Table II. Forecast Levels (cb) 



P (cb) 




p(cb) 




p(cb) 


95.86 


19 


44.81 


37 


20.95 


91.9 


20 


42.96 


38 


20.08 


88.1 


21 


41.18 


39 


19.25 


84.45 


22 


39.48 


40 


18.25 


80.96 


23 


37.85 


41 


17.69 


77.61 


24 


36.28 


42 


16.96 


74.4 


25 


34.78 


43 


16.26 


71.32 


26 


33.34 


44 


15.59 


68.37 


27 


31.96 


45 


14.94 


65.54 


28 


30.64 


46 


14.32 


62.83 


29 


29.37 


47 


13.73 


60.23 


30 


28.16 


48 


13.16 


57.54 


31 


26.99 


49 


12.62 


55.35 


32 


25.88 


50 


12.10 


53.06 


33 


24.81 


51 


11.60 


50.87 


34 


23.78 


52 


11.12 


48.76 


35 


22.80 


53 


10.66 


46.75 


36 


21.85 


54 


10.22 
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